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ABSTRACT 

In this paper we have made an attempt to derive the radial profiles of density 
and temperature of intracluster gas based on the two well-established facts 
at present: the X-ray observed surface brightness of clusters described by the 
standard /3 model and the (NFW) universal density profile as the underlying 
dark matter distribution. We have numerically solved the hydrostatic equation 
by demanding that the volumed-averaged baryon fraction of a cluster should 
asymptotically approach the universal value at its viral radius. We have 
shown that the radial temperature variation derived from these constraints 
differs significantly from the conventional polytropic equation of state: The gas 
temperature profile may show a dramatic increase or decrease with outward 
radius, depending sensitively on the [3 parameter. A large f3 value (typically 
> 0.8) is required in order to ensure that the X-ray temperature makes 
a drop at the virial radius. This indicates that either the NFW profile is 
inappropriate to apply for the overall dark matter distribution of clusters or 
other non-gravitational heating processes may play an important role in the 
dynamical evolution of clusters. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: 
general — X-rays: galaxies 
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1. INTRODUCTION 

One of the major uncertainties in the cosmological apphcations of dynamical properties 
and evolution of galaxy clusters, e.g., the baryon fraction /b and the X-ray luminosity- 
temperature correlation, has probably arisen from the poorly constrained radial variations 
of density and especially, temperature of hot intracluster gas. Indeed, the currently 
available X-ray spectral data of clusters in literature may have suffered various instrumental 
limitations (see, for example, Irwin, Bregman & Evrard 1999), and even conflicting results 
regarding the radial temperature gradients can be obtained by different authors based on 
different analysis techniques (e.g. Markevitch et al. 1998; Irwin et al. 1999; White 2000; 
Irwin & Bregman 2000). This leads to an embarrassing situation that an oversimplified, 
isothermal gas distribution is often adopted in the estimate of the total gravitating mass 
^tot(^) of a cluster enclosed within radius r via the hydrostatic equilibrium hypothesis. 
Here, it is worthy of mentioning another possibility of obtaining the radial distribution of 
intracluster gas in an indirect way, namely, by the inversion of the X-ray and S-Z observed 
surface brightness profiles of clusters (Silk & White 1978; Yoshikawa & Suto 1999). Yet, 
the actual application of this technique to clusters is still limited primarily by the current 
S-Z measurements. 

Several efforts have thus been made towards deriving the radial profiles of density and 
temperature of intracluster gas based on the well-motivated physical mechanisms. Among 
these the shape of the gas density distribution following the universal density profile of the 
dark matter halos of clusters suggested by numerical simulations (Navarro, Frenk & White 
1995; NFW) has received great interest (e.g. Makino, Sasaki & Suto 1998; Ettori & Fabian 
1999; Wu & Xue 2000a; Wu 2000; etc.). While its inner and outer slopes are still under 
debate, the NFW profile may be considered to be a good approximation of the matter 
distribution of dark halos. Alternatively, a number of high-resolution simulations have also 



-4- 



shown that the NFW profile is independent of mass, the initial density fluctuation, and the 
cosmology (Cole & Lacey 1996; Navarro, Frenk & White 1997; Eke, Navarro & Frenk 1998). 
Indeed, within the framework of isothermal and hydrostatic equilibrium there is a striking 
similarity between the gas density distribution predicted by the NFW profile and the 
conventional f3 model revealed by the X-ray observations (Makino et al. 1998). Inclusion 
of the self-gravity of the gas and the polytropic equation of state makes it promising to 
directly link the theoretical predictions with the actual X-ray observations of clusters (Suto, 
Sasaki & Makino 1998). 

Here, we view the problem from a different angle: Given (l)the radial variation of 
the dark matter particles characterized by the NFW profile and (2)the X-ray observed 
surface brightness profiles of clusters, what can we learn about the radial density and 
temperature profiles of the intracluster gas ? Namely, we begin with the two well-established 
'observational' facts today, and drop the assumption about the temperature properties 
of the gas (isothermal or polytropic). Nevertheless, we include the self-gravity of the 
gas. We then solve the hydrostatic equation with the following boundary constraint: The 
volume-averaged baryon fraction fh{r) within the virial radius r^ir should asymptotically 
approach a universal value /b,BBN = ^b/^m defined by the Big Bang Nucleosynthesis, where 
fib and fijii are the average baryon and total mass densities of the Universe in units of the 
critical density pcrit for closure, respectively. The reasons why we impose this condition on 
the solution are as follows: First, it is generally believed that the matter composition of 
clusters averaged over a sufficiently large volume should be representative of the Universe. 
Second, the present X-ray observations have revealed an increase in the baryon fraction of 
clusters with radii and no evidence for an asymptotic tendency towards a universal value 
at large radii (White & Fabian 1995; David 1997; White, Jones & Forman 1997; Ettori & 
Fabian 1999; Jones & Forman 1999; Markevitch et al. 1999). This conflict has probably 
arisen from the oversimpliflcation of the currently adopted temperature model for the 
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intracluster gas. Setting a priori the universality of /b(^vir) can instead allow one to work 
out the constraints on the radial variations of density and temperature of the gas. Third, 
the boundary condition will be required mathematically in order to solve the hydrostatic 
equation, while little has been known about the gas properties within the central cores of 
clusters where the local dynamical activities may play a dominant role. 

In a word, if the distribution of dark matter can be reasonably approximated by 
the NFW profile, we will be able to derive the radial density and temperature profiles 
of intracluster gas based on the well measured X-ray surface brightness of clusters with 
the constraint that the baryon fraction within the virial radius is universal. The resulting 
density and temperature profiles can be directly compared with the X-ray spectroscopic 
measurements and the S-Z observations of clusters, which will provide a critical test for 
the validity of the NFW profile and the hydrostatic equilibrium hypothesis. Meanwhile, 
such a study will be helpful for resolving the conflict regarding the different temperature 
measurements, i.e., whether or not the X-ray emitting gas can be regarded as isothermal in 
clusters. Eventually, we hope that this research will be of significance for our understanding 
of the fundamental physical processes in the dynamical evolution of clusters including 
non-gravitational heating. 

2. THE MODEL 

We assume that the dark halo of a cluster is described by the NFW profile 

Pdm(^) = / , w-i , — Tyl-' (-*■) 

(r/rs)(l + r/r^Y 

where ps and r^ are the characteristic density and length, respectively. Another two 
equivalent parameters frequently used in literature are 5c = Ps/Pait and c = r^ir/rs- The 



-6- 

virial radius r^i^ is defined by 

47r 

MY)Mij\ir) = — r^ij,AcPcrit, (2) 

where Ac represents the overdensity of dark matter with respect to the average background 
value pcrit, for which we will take Ac = 200 in the following computation. Our main 
conclusion is unaffected by this choice. The total mass of the dark matter halo enclosed 
within radius r is 



Mdm(^) = 4:71 p J- . 



In 1 



(3) 



ngas(r)TV4(r) = n,^X^' 1 + ^ , (5) 



Tg / r + r, 

On the other hand, the X-ray surface brightness profile S'x(r) can be approximated by the 
conventional (3 model: 

/ ^2 \ -3/3+1/2 

^x(r) = So[i + ^J , (4) 

in which Tc denotes the core radius. For the thermal bremsstrahlung emission, the above 
form of Sx{r) indicates (Cowie, Henriksen & Mushotzky 1987) 

,2 \ -3/3/2 

where Ugas and T are the gas number density and temperature, respectively. Equation (5) 
identifies the f3 density profile in the case of isothermality. As a result, the total mass in gas 
within r is simply 

Mg^{r) = 47r/impngas,o J [yj ( "^ + ~ J ^^^^' ^^^ 

in which n = 0.585 is the average molecular weight. The volume-averaged baryon fraction 
/b(r) within r is defined as 

^'^ ^ Mgas(r) + MDM(r)- ^' 

Here we have neglected the contribution of stellar mass to /b(^)- We further assume that 
the gas is in hydrostatic equilibrium with the underlying gravitational potential of the 
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cluster produced by Mtot(^) = Mgas(^) + Mdm(^), i-e., 

G'[Mgas(r) + MDM(r)] _ 1 d{n^^kT) 

r^ /xmpngas dr 

Using temperature and baryon fraction as the two variables and re-organizing the above 
equations yield 

df A(3xf Aa, 



dx x"^ + a^ 3a;2 



ln(l + x)- 



l + x 



l-/b' 



(9) 



dx in(l+x)-^ ' ^ ' 

in which the scaled quantities are: T = T/Tq, x = r/vs, a = Vc/rs, b = /impngas,o/Ps, and 
ao = AttG firripPsr'^ / kTo. Our task is thus reduced to finding the solutions of equations (9) 
and (10) with the following boundary conditions 

f (0) = 1, (11) 

and 

fb{c) = /b,BBN; (12) 

f^ U=c = 0. (13) 

ax 

This last constraint demands that the baryon fraction should asymptotically match the 
universal value of /b,BBN at the virial radius r^h or equivalently c. 



3. RESULTS 

3.1. General properties 

There are a total of five parameters involved in equations (9) and (10): j3, ao, a, b 
and 6c or c. Additionally, the universal baryon fraction /b bbn should also be regarded as 



another unknown parameter although fib has been fixed by the Big Bang Nucleosynthesis. 
However, these parameters are not independent, as we will discuss in some detail later. 
Mathematically, the number of these free parameters can be reduced by two with the 
restriction of the two first-order differential equations, equations (9) and (10), and the 
associated boundary conditions. Namely, given four of the six parameters (/3, aQ, a, b, 6c, 
/b,BBN)) we will be able to find the other two from equations (9) and (10). Technically, 
we perform iteratively the numerical searches for the solutions of equations (9) and (10) 
over a space of two parameters until the boundary conditions equations (11)-(13) are 
satisfied. Using the resultant temperature T, we will then derive the scaled gas density 
figas = '^gas(^)/%as,o accordiug to equation (5): 

/ ^2 \ -3/3/2 

n,as = T-V4 (^1 + _j . (14) 

We demonstrate in Figure 1 a set of typical solutions for a = 0.3, 6c = 10"^ and /b,BBN = 0.1 
but with different choices of (3: 2/3, 0.75, 0.85 and 1. The other two parameters, b and aQ, 
can be uniquely determined in each case during our numerical searches for the solutions, 
which read (6, ao) = (0.34, 5.73), (0.42, 6.22), (0.65, 6.61) and (1.34, 6.96) for /3=2/3, 
0.75, 0.85 and 1, respectively. The halo concentration for our example is found to be 
c = r^ir-frs = 5.32. We have also displayed in Figure 1 the input 5'x(r) and our derived 
radial profiles of gas density, Ugas- 

EDITOR: PLACE FIGURE ?? HERE. 

It is not surprising from Figure 1 that the shape of the resulting radial profiles of gas 
density, n^^, is very similar to the conventional expectation, i.e., the (3 model. Alternatively, 
the baryon fraction fh{r) shows an increase with radius roughly within the core regions, 
and then asymptotically approaches the background value /b,BBN at virial radius c (for a 
small (3) or reaches a maximum before it matches the background value (for a large f]). 
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However, the temperature profile T{r) derived from equations (9) and (10) along with the 
restrictions equations (11)-(13) exhibits a remarkable difference in the radial variation 
among different (3 clusters, from a rapid increase at large radii for a small /? to a sharp drop 
near the virial radius for a large f3, separated roughly at /3 ~ 0.8. This property is mainly 
due to our restriction on the baryon fraction at a; = c (equations [12] and [13]): /b = /b,BBN 
and dfb/dx = 0. Recall that a similar temperature behavior is required in order for fi, to 
asymptotically approach the universal value /b,BBN at large radii within the framework of 
the f3 model and the hydrostatic equilibrium for intracluster gas, in which the separation of 
the two different variations of the asymptotic temperature at large radii occurs ai (3 = 2/3 
(Wu & Xue 2000b). 

Similarly, we can numerically obtain the eigenvalues (ao, b) for different choices of {f3, 
0-1 ^ci /b,BBN)- In Figure 2 we display the dependence of the resultant eigenvalues (ao? b) 
upon each of the four parameters. The most interesting result is the cut-off of j3 parameter 
set by ao = 0, below which, the unphysical solutions will occur. This indicates that the 
observed X-ray surface brightness profile S^{r) described by the (3 model should not have 
too small values of (3 if our working scenario holds true. In order to demonstrate how the 
cut-off value, /Scut, varies with clusters, we perform a search for the solutions with ao = 0. 
In this case, the eigenvalue h can be uniquely fixed by our numerical solutions. Meanwhile, 
our numerical search shows that /b,BBN produces a minor effect on /3cut- Consequently, /Scut 
depends mainly on two parameters: 5c and a (Figure 3). It turns out that for a typical 
cluster with 5c = 10^-10^ and a = 0.1-1, the (3 parameter should be greater than 0.5, and 
/Scut Hiay exceed 0.6 if 5c > 10^. These constraints on j3 parameter seem to be in good 
agreement with the X-ray observations that the majority of the X-ray clusters have X-ray 
surface brightness profiles that satisfy j3 > 0.5. 

EDITOR: PLACE FIGURE ?? HERE. 
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Essentially, the parameter ao depends very weakly on Sc, /b,BBN and P (for P > /5cut), 
and its typical values for a = 0.15 and 0.3 are around 10 and 6.5, respectively. Our estimates 
can be compared with the previous studies that assumed an isothermal gas distribution 
tracing the gravitational potential of the NFW profile but dropping the restrictions of 
equations (12)-(13): For rich clusters the variation of ao is limited to a very narrow range 
from 6 to 25 with an average value of ~ 10 (Ettori & Fabian 1999; Wu & Xue 2000a; Wu 
2000). Actually, ao can be approximately regarded as a constant for a given a (see Figure 
2). 

3.2. Asymptotic temperature at large radii 

We learn from Figure 1 that there exists a critical value of /3crit such that T{c) = 1. 
Namely, the radial temperature profile close to x = c will show a tendency to increase 
(decrease) with radius when P < /3crit {P > Pait)- The existence of such a /3crit can be 
understood qualitatively by considering the following limiting situation where the virial 
radius c is much greater than the gas core radius a. The universal profile admits an x~^ 
density at large x. If the distribution of /b,BBN is to be asymptotically flat, the gas density 
Ugas must also obey x^^ at large x. Given that the P model requires equation (14) to hold, 
the existence of a Pent thus naturally follows. For this limiting situation, /Scrit = 1 but when 
c no longer signiflcantly exceeds a, the value of /5crit is found to vary somewhat as described 
below. 

While the asymptotic high-temperature and low-temperature are both allowed 
mathematically, there should be one temperature proflle at large cluster radii. Indeed, it is 
very unlikely that the gas temperature would increase monotonically with outward radius. 
At least, there is no known mechanism that can produce a signiflcant rise of gas temperature 
at outermost radii. Note that this does not exclude a mild increasing temperature with 
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radius over a certain region of clusters (e.g. Irwin & Bregman 2000). As a result, the 
condition T{c) = 1 sets another low limit on the value of f3, which guarantees that the gas 
temperature would not rise at virial radius when j3 > /dent- 
in order to determine /5crit, we once again perform the numerical searches for the 
solutions of equations (9) and (10) with the additional constraint T{c) = 1. This last point 
reduces the free parameters to three: /b,BBN, ^c and a. In Figure 3 we plot the resulting /5crit 
against 6c for two different choices of a (a = 0.1 and a = 1) but a fixed universal baryon 
fraction of /b,BBN = 0.1. It appears that although /Scrit shows a radial variation similar to 
that of /3cut, the value of /?crit is apparently larger than the corresponding result for /3cut, 
indicating that the observed X-ray surface brightness profiles of clusters should possess an 
even larger (3 parameter than jScut to ensure the natural drop of gas temperature at the 
edges of clusters. For our examples shown in Figure 3, the (3 parameters are required to 
exceed 0.7 and 0.8 for 6c = 10^ and 10'^, respectively. This may constitute a critical test for 
the working hypotheses established in the above section. The currently available data from 
X-ray measurements seem to marginally agree with these limits even if the effect of cooling 
fiows (e.g. Vikhlinin et al. 1999), and the unknown parameters (a, 6c and /b,BBN) are taken 
into account. We will explore in detail the possible reasons in the discussion section. 

EDITOR: PLACE FIGURE ?? HERE. 



3.3. Self-consistent solution 

In an optically thin plasma emission model, the central gas density ngas,o is related to 
the central surface brightness 5'x(0) through 

^^^'°" ^7Ia(To)r(3/?-l/2) r^ ' ^''^ 
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where q;(To) is the cooling function 

'g is the average Gaunt factor and JI = yU^(l + X)/2 with X = 0.768 being the primordial 
hydrogen mass fraction. The X-ray measurement of the surface brightness profile of a 
cluster provides us with the following quantities, S'x(O), /J and Vc- On the other hand, the 
central electron temperature Tq, can be obtained relatively easily by the X-ray spectroscopic 
measurements, as compared with the outermost temperature, e.g. T{x ~ c). Given S'x(O), 
j3, Tc and Tq, we will be able to derive the central gas density ngas,o in terms of equations 
(15) and (16). Moreover, since the purpose of the present paper is not to estimate the 
baryon fraction of clusters, we assume that the universal value of /b,BBN is known. As a 
result, there are only two parameters, p^ (or equivalently 5c) and r^ in the NFW profile, 
that remain to be determined. We thus search for the solutions of equations (9) and (10) 
by iterations with a set of (5c, a) until a self-consistent solution is achieved, which uniquely 
determines the parameters 5c and r^. in the NFW profile. Recall that in the isothermal case, 
the X-ray surface brightness profile predicted by the NFW profile is directly applicable to 
the fitting of the observed >S'x(r), giving rise to the two parameters, 5c and r^ (Makino et 
al. 1998; Ettori & Fabian 1999; Wu & Xue 2000a,c; Wu 2000). Indeed, it appears that the 
procedure of finding the two parameters in the NFW profile now becomes very complicated 
if the temperature profile of a cluster is unknown. 

In the present paper we do not intend to determine the parameters (ps, r^) by applying 
the current method to an ensemble of X-ray clusters. Instead, we present the result for 
a typical cluster with ngas,o = 5 x 10~^ cm~^, r^ = 0.25 Mpc {Hq = 50 km s~^ Mpc^-*^) 
and To = 7 keV but different [3 parameters. We fix the universal baryon fraction to 
be /b,BBN = 0.1. The resulting ps and rg, along with the concentration parameter c (or 
equivalently, the virial radius) and the corresponding temperature T{c) at a; = c, are listed 
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in Table 1. Apart from the significant differences in the outermost temperatures at virial 
radii for different {3 parameters, the values of the remaining parameters, rg, 5^ and c, all 
fall inside the limits for typical rich clusters with an isothermal gas distribution (Ettori 
& Fabian 1999; Wu & Xue 2000a,c; Wu 2000). This implies, on the other hand, that the 
determination of the two parameters in the NFW profile from the X-ray observations of 
intracluster gas is probably insensitive to the radial temperature variations. However, the 
temperature rise or drop near the edge of a cluster depends critically on the (3 parameter. 

4. DISCUSSION AND CONCLUSIONS 

The current X-ray measurements of the radial temperature profiles of intracluster 
gas especially at outer radii have some ambiguity, which is probably responsible for the 
major uncertainties in the present determination of the gravitating masses of clusters. 
Additionally, the large dispersion of the derived baryon fraction among different clusters 
may be also associated with the poorly constrained cluster temperatures. Consequently, the 
cosmological application of the dynamical properties of clusters are significantly affected. 
In the absence of a reliable constraint on the radial temperature variation from X-ray 
spectral measurements, we have made an attempt in the present paper to derive the 
temperature profile T{r) based on the well-established facts from the X-ray observations 
and numerical simulations, along with some plausible boundary conditions. Specifically, 
we use the NFW profile as the dark matter distribution udmI^) of clusters, and assume 
that intracluster gas ngas(^) is in hydrostatic equilibrium with the underlying gravitational 
potential dominated by nDM(^) and n^g^^{r). We adopt the X-ray observed surface brightness 
profile S^{r) described by the conventional (3 model to set up a link between gas density 
ngas{r) and temperature T{r). We have then numerically solved the hydrostatic equation 
by demanding that the baryon fraction /b(r) approaches asymptotically the universal value 
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/b,BBN = ^b/^M at the viral radius. As a result, we can simultaneously determine the 
radial profiles of gas density, temperature and baryon fraction in clusters. In particular, 
using the X-ray data of /3, Vc, >S'x(0) and Tq, we are able to find the two free parameters in 
the NFW profile, with which the dark matter profile can be completely fixed. 

The temperature profile constructed jointly from the NFW profile and the X-ray 
surface brightness measurement appears to be significantly different from the conventional 
speculation. Apparently, our temperature profiles (see Figure 1) cannot be simply described 
by the widely adopted polytropic equation of state, T oc n'^^^ where 7 is the polytropic 
index, which is refiected by the striking difference between the core radius of gas density 
and that of temperature. It seems that the shape of temperature profiles in the outer 
regions of clusters can be dramatically different, depending sensitively on the (3 parameter. 
This point has mainly arisen from the restriction that the volume-averaged baryon fraction 
should asymptotically approach the universal value at the virial radius. The recent study by 
Wu & Xue (2000b) based on the (3 model and the polytropic gas distribution has essentially 
reached a similar result, i.e. a mild increase in temperature with radius should be required 
for clusters with 1/3 < j3 < 2/3 in order to ensure the universal constancy of the cluster 
baryon fraction at large radii. Although this possibility has been reported for some clusters 
by Irwin & Bregman (2000) and shown in Figure 1 for small (3 (more clearly Figure 4 for 
f3 = 2/3), naturally we would expect a drop in temperature at the outermost radii. This can 
be achieved if the /3 parameters become larger (typically (3 > 0.8). While numerical studies 
of cluster formation and evolution often yield a similar result of /3 ~ 0.7-1, observationally, 
it is unlikely that one can significantly raise the f3 value to j3 > 0.8 from the fitting of the 
X-ray surface brightness profiles of clusters. Exclusion of the central cooling regions or 
employment of the double f3 models in the fitting of the observed Sx{r) can only lead to a 
moderate increase of the f3 parameter (e.g. Jones & Forman 1984; Vikhlinin et al. 1999; 
Xue & Wu 2000a). 



-15- 

In addition to the X-ray spectroscopic measurements, the SZ effect may also help to 
make a distinction between the temperature rise or fall models at large radii close to r^ij.. 
The SZ effect depends on the gas density and temperature via AT/Tcbr oc / rigssTdi, where 
TcBR represents the temperature of the cosmic background radiation, and the integral is 
performed along the line of sight i. In Figure 4 we demonstrate the thermal SZ effect 
produced by two clusters with f3 = 2/3 and 1, respectively, adopting the same parameters 
as in Figure 1. If the gas temperature shows a dramatic increase with radius, the shape 
of the SZ profile will be strongly distorted. On the contrary, it turns out to be hard to 
distinguish the temperature decreasing model from the isothermal one simply based on the 
SZ observations. 

EDITOR: PLACE FIGURE ?? HERE. 

Although we have fulfilled our task that the radial profiles of density and temperature 
of intracluster gas are derived from the well-established facts, the existing X-ray data 
and future X-ray observations can provide a critical test of these predictions. Any 
conflict between observations and our predictions would imply that at least one of the 
following working hypotheses should be abandoned: (l)the NEW profile as the dark matter 
distribution of clusters; (2)the hydrostatic equilibrium for intracluster gas; (3)the pure 
thermal bremsstrahlung as the observed X-ray emission; and (4)the total baryon fraction 
of clusters as a reliable indicator of matter mixture of the Universe. In fact, all these 
hypotheses have been challenged in recent years. 

It has been noticed that the central cusp in the NEW profile disagrees with the soft 
inner mass distributions of galaxies and clusters, which has led Spergel & Steinhardt (2000) 
to propose the possibility that the cold dark matter particles are self-interacting. Indeed, 
the X-ray surface brightness of clusters predicted from the NEW profile has a much smaller 
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core radius than the one seen from X-ray observations (Makino et al. 1998). In fact, the 
cusped NFW profile yields a higher central gas temperature, which can be clearly seen 
from the temperature profiles plotted linearly in Figure 4, while the X-ray spectroscopic 
measurements often reveal a cold gas component inside the core regions. As a result, a 
modification to the NFW profile should be properly made, such as the empirical density 
profile with a definite core radius Pdm = Po(l + x)^^{l + x^)^^ (Burkert 1995), which may 
both provide a good fit to the dark matter distribution and predict the right shape of 
the X-ray surface brightness distribution of clusters (Wu & Xue 2000c). However, fitting 
the predicted X-ray surface brightness profiles by such a revised dark halo profile to an 
ensemble of the X-ray observed surface brightness profiles of clusters, Wu & Xue (2000c) 
have shown that the core radii of the dark matter halos reach the X-ray core sizes of ~ 0.2 
Mpc. The recent numerical simulations under some simple models for self-interacting dark 
matter particles have essentially arrived at the similar conclusion (Yoshida et al. 2000). 
Such large dark matter cores apparently exceed the ones required by modeling giant arcs 
seen in the central regions of clusters (Hammer 1991; Wu & Hammer 1993; Grossman & 
Saha 1994; Tyson et al. 1998). In a sense, any attempts at significantly modifying the 
NFW profile may require novel physical mechanisms. 

Alternatively, it has been argued in recent years whether the intracluster gas is in 
perfect hydrostatic equilibrium with the underlying gravitational potential, based mainly 
on the presence of substructures and complex temperature patterns in clusters revealed 
by optical/X-ray observations (e.g. Henriksen & White 1996; Markevitch 1996). Yet, 
a statistical comparison of cluster masses derived from optical/X-ray observations and 
the NFW profile, under the assumption of hydrostatic equilibrium, and the strong/weak 
gravitational lensing measurements shows excellent agreement among different cluster mass 
estimates on scales greater than the X-ray core radii (Allen 1998; Wu et al. 1998; Wu 2000; 
references therein). This suggests that the hydrostatic equilibrium may break down only 
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on small scales where the local dynamical activities make a non-negligible contribution to 
cluster mass estimates, and hydrostatic equilibrium can be reasonably applicable to the 
overall distribution of intracluster gas. 

Our emphasis is thus laid on the effect of non-gravitational heating processes in the 
dynamical evolution of clusters. It has been known in the past years that the observationally 
determined X-ray luminosity (Lx) and temperature relation, L^ oc T^, deviates significantly 
from the prediction, Lx oc T^, under the standard scenario that assumes the hydrostatic 
equilibrium and the thermal bremsstrahlung emission for intracluster gas (e.g. David et 
al. 1993; Markevitch 1998). While one cannot definitely exclude the possibility that the 
baryon fraction of clusters may vary with temperature (David et al. 1993), this conflict 
is strongly suggestive of the importance of the non-gravitational heating processes in the 
cosmic evolution of clusters. In fact, supernovae in cluster galaxies would substantially 
inject energy into intracluster medium through galactic winds in the early phase of cluster 
formation, and the intracluster gas can be heated to temperatures even exceeding that 
produced by gravitational heating alone (White 1991; David, Forman & Jones 1991). 
This scenario may allow one to reproduce naturally the observed Ly^-T relation L^ oc T^ 
(Ponman, Cannon & Navarro 1999; Wu, Fabian & Nulsen 1999; Loewenstein 2000). In 
particular, the energy input from supernova-driven protogalactic winds will permit the 
intracluster gas to extend out to larger radii, giving rise to a shallower X-ray surface 
brightness profile (David et al. 1990; Ponman et al. 1999; Lloyd-Davies, Ponman & Cannon 
2000). The discrepancy between the large (3 parameters {13 > 0.8) required to maintain the 
decreasing temperature profile near virial radius discussed in the present paper and the 
relatively small value of /5 ~ 0.7 found from the X-ray observed surface brightness profiles of 
clusters is probably due to our negligence of the effect of other heating processes. A reliable 
construction of the radial density and temperature profiles of intracluster gas should also 
allow the non-gravitational heating to be included. 
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Taking these arguments as a whole, we feel that the exercise made in this paper, 
combined with the X-ray observations of clusters, will be helpful for our understanding 
of the fundamental physical process that dominates the dynamical evolution of clusters. 
Applications of the current method to real X-ray clusters will be presented in subsequent 
work (Xue & Wu 2000b). 
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of the National Taiwan University, where part of this research was carried out. This work 
was supported by the National Science Foundation of China, under Grant 1972531, and the 
National Science Council of Taiwan, under Grant NSC89-2112-M008-037. 
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Table 1: 6c and rg for a typical cluster'' 



p 


r, (Mpc) 


^^c 


c 


Tic) 


0.70 


0.51 


7051 


4.54 


4.66 


0.75 


0.49 


7314 


4.63 


2.56 


0.80 


0.41 


8068 


4.84 


2.00 


0.85 


0.39 


8600 


4.98 


0.77 


1.00 


0.29 


11000 


5.54 


0.06 



-3 ovvi-3 



To = 7 keV, re = 0.25 Mpc and Ugas.o = 5 x 10"^ cm 
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Fig. 1. — The derived radial profiles of gas density (ngas/%as,o), temperature (T/Tq) and 
baryon fraction (/b) for /b,BBN = 0.1, a = Vc/rg = 0.3 and 6c = Ps/Pcrit = 10^ but with four 
different choices of j3. Also shown are the input X-ray surface brightness profiles of clusters 
described by the conventional f3 model. The dotted line denotes the virial radius. 
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Fig. 2. — Dependence of q;o(= 4:7iG firripPsr^ / kT^) (dotted lines) and b{= /impngas,o/Ps) (solid 
lines) on (1)/? (5^ = 10^ a = 0.3 and /b,BBN = 0.1), (2)/b,BBN {P = 0.75, 5c = 10^ and 
a = 0.3), (3)4 {p = 0.75, a = 0.3 and /b,BBN = 0.1), and (4)a (/5 = 0.75, 4 = 10^ and 

/b,BBN = 0.1). 
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Fig. 3. — Lower limits on j3 parameter against the characteristic density 6c for two choices 
of a: (l)a = 0.1 (sohd hne and open squares) and (2)a = 1 (dashed hne and open circles). 
Solid and dashed lines: the limits /3cut from the unphysical solutions of equations (9) and 
(10) as a result of ao = 0; Open squares and circles: the limits /3crit set by the requirement 
of T(c) /To = 1. 
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Fig. 4. — A comparison of the SZ profiles predicted by tlie isothermal gas distribution 
(dotted lines) and our derived radial distribution of intracluster gas from equations (9) and 
(10) (solid lines). Upper and lower panels correspond to the situations of /9 = 2/3 and /9 = 1 
in Figure 1, respectively. The temperature profiles are shown by dashed lines. For clarity, 
the vertical scales for /5 = 2/3 (upper panel) have been multiplied by a factor of 10. 



